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ABSTRACT 

We consider large-scale collective motion of flat edge-on spiral galaxies from the Revised 
Flat Galaxy Catalogue (RFGC) taking into account the curvature of space-time in the Lo- 
cal Universe at the scale lOO/i -1 Mpc. We analyse how the relativistic model of collective 
motion should be modified to provide the best possible values of parameters, the effects that 
impact these parameters and ways to mitigate them. Evolution of galactic diameters, selec- 
tion effects, and difference between isophotal and angular diameter distances are inadequate 
to explain this impact. At the same time, measurement error in H I line widths and angular 
diameters can easily provide such an impact. This is illustrated in a toy model, which allows 
analytical consideration, and then in the full model using Monte Carlo simulations. The re- 
sulting velocity field is very close to that provided by the non-relativistic model of motion. 
The obtained bulk flow velocity is consistent with ACDM cosmology. 
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1 INTRODUCTION 

At present time the Universe is essentially inhomogeneous on the 
scales of about 10-100 Mpc. The development of initial fluctua- 
tions led to an observable large-scale structure. The regions with 
increased matter density provide an additional attraction of sur- 
rounding galaxies. The regions with decreased density, e.g. voids, 
also make an input to the collective large-scale motion of galaxies 
on the background of Hubble expansion. Investigation of such mo- 
tion on one side allows to map the matter density, including dark 
matter, in the Local Universe, and on the other side its parame- 
ters are linked with cosmological parameters. All of this makes the 
study of collective galaxy motions important. 

In recent years a number of articles was published claim- 
ing that typical velocities of large-scale collective motions are 
inconsistent with th e standard ACDM model. For example, 
IWatkins et ail d2009h obtained the value 407 ± 81kms _1 at 
the scale 100/i -1 Mpc, whereas the ACDM model gives about 
250 km s -1 . However, our estima tion of 210 ± 86kms~ at th e 
same scale, obtained in the article (Parnovsky & Parnowski 2010h . 
is consistent with the ACDM predictions. Additionally, in the same 
article we obtained from the peculiar velocities the constraints 
on the cosmological parameters Q m and erg and their combina- 
tions, which match the other more precise constraints like baryonic 
acoustic oscillatio ns or WMAP observations. 

In the article dParnovskv & Parnowski 2010) we used a sam- 
ple of RFGC galaxies with measured redshifts and H I line widths. 
The Revised Flat Galaxy Catalogue (RFGC) fearachentse v et all 
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1 19991) and its previous v ersion Flat Galaxy Catalogue (FGC) 
( Karachent sev et all Il993h contain the information about Right 
Ascension and Declination for the epochs J2000.0 and B 1950.0, 
galactic longitude and latitude, major and minor blue and red di- 
ameters in arcminutes in the POSS-I diameter system, morpho- 
logical type of the spiral galaxies according to the Hubble clas- 
sification, index of the mean surface brightness and some other 
parameters, which are not used in this article. The RFGC con- 
tains data about 4236 flat edge-on spiral galaxies, almost uni- 
formly covering the celestial sphere and satisfying the conditions 
db/bb 7 and a& > 0.'6. Here and bb are the major and mi- 
nor axial diameters in the 025 system. The original goal of this 
catalogue was to estimate the distance to galaxies according to 
the Tully-Fisher relation in the "Hi line width - linear diame- 
ter" version without using their redshifts. The data about the red- 
shifts and H I line widths or gas rotation velocities V ro t were taken 
from different sources. There were a nu mber of gradually increas- 



ing samples of galaxies with such data 1 Karachentse v et al1l2000l : 
IPamovskv et alj|200ll ; IParnovskv & Tugay||2004h . The latest ver- 



sion of thi s sample including 1623 galaxi es was compiled and de- 
scribed bv lParnovskv & ParnowskilbOld) . A list of peculiar veloc- 
ities based up on this list in the non-rela t ivistic model of motion was 
assembled bv lParnovskv & Parno wski (200j|). 



In this article we use the same sa mple, but with differ- 



ent model o f collective motion of galaxies I lKudrva & Alexandrovl 
I2002L |2004|) . based upon the general theory of relativity (GTR). 
This model was applied earlier to the p revious version of the sam- 
ple bv lParnovskv & Gavdamakal d2004h : however, the present arti- 
cle offers a much more in-depth analysis. 
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2 DESCRIPTION OF MODELS OF COLLECTIVE 
MOTION OF GALAXIES 

We assume that there is a three-dimensional velocity field of col- 
lective galactic motions on the background of Hubble expansion. 
We consider the galaxies in the sample not as massive objects, but 
rather as test particles, whose peculiar velocities are indicators of 
the large-scale velocity field. Using the multipole decomposition 
of large-scale velocity field up to quadratic terms in distance and 
switching to the radial component we obtain the expression for ra- 
dial velocity of an individual galaxy. The actually measured radial 
velocities differ from those predicted by this model due to devia- 
tions from the statistical Tully-Fisher relation, influence of motions 
with scales much less than the sample depth (fall towards nearby 
attractors) and random errors. Treating these errors as stochastic 
we can use the least squares method to calculate the parameters of 
the model. 



where W is a corrected H I line width in km s" 1 measured at 50 per 
cent of the maximum, a is a corrected major galaxies' angular di- 
ameter in arcminutes on red POSS and ESO/SERC reproductions, 
U is a ratio of major galaxies' angular diameters on red and blue 
reproductions, T is a morphological type indicator (T = I t — 5.35, 
where I t is a Hubble type; It — 5 corresponds to type Sc), and 
B is a surface brightness indicator (B — Isb — 2, where Isb is 
a surface brightness index from RFGC; brightness decreases from 
I to IV). Note that the statistical significance of eac h term in eq . 
(|6) is greater than 99 per cent according to the F-test jFisherllT954l ; 
lHudson|[l964l) . 

Thus, the DQO-model contains 24 parameters, namely 3 com- 
ponents of the vector D, 6 coefficients d, 5 parameters g; of the 
tensor Q, 3 components pi of the vector P and 7 parameters o, of 
the tensor O. By rejecting V oct we get a simpler DQ-model with 14 
components. Further rejecting V qua leads to the simplest D-model 
with 9 components. 



2.1 Non-relativistic model of collective motion 

Before discussing the relativistic model of collective motion, let 
us briefly recall the non-re lativistic m odels, introduced in the arti - 
cle I Parnovsky et al.l200ll) and used bv lParnovskv & Tugavl fe004h 
and Parnovsky & ParnowskJ d201(j|) . We will start from the more 
complex of them, namely the DQO-model. 



V = R + V dlp + V qua + V oct + SV. 



(1) 



Here Visa radial velocity of the galaxy in the CMB isotropy frame, 
R — Hr is the Hubble velocity, r is the distance to the galaxy, 
SV is a random error, V dip , V qua and V oct are the dipole (D), 
quadrupole (Q) and octopole (O) components of the non-Hubble 
cosmic flow. They are given by the following expressions: 

V dip = Dim, (2) 

V qua = RQ ik rHn k = R ( qi (n( - n 2 3 ) + q 2 (n\ - n%) 

y OC t = R 2( 0ikintnkni + p i7l -) = R 2 (p^. 

+ oi(3nm| - nf) + o 2 (3nm 3 — n\) 

+ o 3 (3n 2 n? - n|) + o 4 (3n 2 n§ - n%) (4) 

+ o 5 (3n 3 n? - n 3 ) + o 6 (3n 3 nl - n 3 ) 

+ o 7 nin 2 n 3 ) . 

Here we use the Einstein rule - summation by repeated indices; rn 
are Cartesian components of the unit vector n towards the galaxy, 
connected with the galactic coordinates / and b in the following 
way: 



n\ — n z — sin b, 
n 2 — n x = cos / cos b 
ri3 



(5) 



n y — sin / cos b. 



The dipole component (bulk motion) is described by the vector D. 
The quadrupole component is described by the symmetrical trace- 
less tensor Q. It has 5 independent parameters qi. The octopole 
component can be described by one rank 3 tensor, but we divide 
it into a trace characterized by vector P and a tensor O, which is 
antisymmetrical with respect to each pair of indices. The latter has 
7 independent parameters Oi . 

Hubble velocity is determined from the generalized Tully- 
Fisher relation in the "angular diameter -Hi line width" version 
by the following formula 

, w 



R=(d+ C 2 B + C Z BT + dU)^ 



(6) 



2.2 Relativistic model of collective motion 

The existence of multipole components is due to density inho- 
mogeneities in the Universe. For homogeneous isotropic Universe 
instead of equation QJ we would have the simple Hubble law 
V = Hr. Hubble expansion is due to cosmology and thus is most 
adequately described in the framework of GTR. Such a descrip- 
tion raises a problem of distance measurement. The main types of 
distances used are the redshift distance, photometric distance and 
angular diameter distance. They are related to each other through 
formulae containing cosmological parameters. In our case the nat- 
ural choice is the angular diameter distance, since we determine 
distances by the generalized Tully-Fisher relation using angular di- 
ameters. Thus, when speaking about distance r in relativistic mod- 
els, we will mean angular diameter distance. 

For the homogeneous isotropic cosmological models the de- 
pendence of the velocity V = cz on R = Hr in the next order in 
R has the form 



V = R + jR 2 



(7) 



The coefficient 7 is expressed through the deceleration parameter q 
by 

3 + g 



7 : 



2c 



(8) 



For different types of distance another expression for 7 should be 
used. For the standard ACDM cosmology we have 

q = _S _ fi A = _o.61, (9) 



where Q m and S7a are relative densities of matter, including dark 
matter, and dark energy respectively. Numerical estimations are 
based on the results of 7-year WMAP observations dKomatsu et al.l 
2010). Introducing l[9} into l[8} we obtain 

3.98- 10" 6 skm _1 . 



(10) 



Now let us consider a relativistic model of galaxy 
motion in inhomogeneous space -time. It was developed by 
iKudrva & Alexandre^ j2002Ll2004l) . Instead of {D we use the equa- 
tion 



V rel =R+ V d ' p + V qua + V° ct + jR 2 + SV. 



(11) 



Here V re l is still given in the CMB frame, the expressions for the 
dipole l[2j and qudrupole l[3j components remain the same, and the 
octopole component assumes the form 
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V°2i = R 2 {Pirn + Oijkrurijnk + Sijmrij 



(12) 



Here S is a symmetric traceless tensor, characterized by 5 indepen- 
dent parameters sc 

StjTurij = si(nf - 7T-1) + s 2 (nj - n%) 

+ S3mri2 + S4,n\ri3 + s^^ns. 



(13) 



The parameter 7 is related to the deceleration parameter by 



7 = 



3 + g 
2c 



— QikQ lk - 
3c 



(14) 



This equation reduces to equation {8} in absence of quadrupole 
component, for example, for homogeneous models. Similarly to the 
non-relativistic case, we can reduce the relativistic DQO-model to 
DQ- and D-models b y rejecting highest-order mult i poles. 

In the papers jKudrva & Alexandre^ I2002L |2004) it was 
shown that the relation 



a 



MS 



V?V S = %r {-2cS al + 6Q Q7 - 3Q%Q ej 
- (V a Vj - g ai Qe V Q tv )) 



(15) 



must hold. Here Greek indices denote four-dimensional compo- 
nents, g a p is a metrical tensor, and C a p 1 s is a Weyl tensor. The 
spatial parts of four-dimensional tensors Qap and S a p coincide 
with three-dimensional tensors Q and S, and their temporal and 
mixed parts are much smaller. The four-dimensional velocity vec- 
tor V a can be assumed equal to ((fa, , 0, 0, 0). 

Now all that remains is to supply the relativistic model with an 
expression for the angular diameter distance. Since we consider the 
terms proportional to R 2 separately, we should remove the terms 
quadratic in distance from the generalised Tully-Fisher relation {6j: 

W 1 

R= (Ci + C2B + C3BT + C4U)— +C 5 -. (16) 

a a 

Note that all the remaining terms in this equation are inverse pro- 
portional to the angular diameter a. 

Using the data about radial velocities, H I line widths, angu- 
lar diameters, morphological types, surface brightness indices and 
axial ratio of galaxies it is possible to obtain the value and er- 
rors for all 29 parameters of the relativi stic model. For the previ- 
ous ve rsion of the sample it was done bv lParnovsky & Gavdamakal 
d2004) . The statistical weights of all galaxies are assumed to be 
the same. Besides the whole sample we also use subsamples lim- 
ited in depth R < Rmax- They are defined as follows: they con- 
tain all galaxies, which have Hubble velocity less than R max in 
the non-relativistic D-model. The results of processing of the sub- 
samples with Rmax = 8000 km s -1 and Rmax = 10000 km s _1 
are presented in Table [T] It contains information about the stan- 
dard deviation a, coefficients of the Tully-Fisher relation, apex 
coordinates and modulus of the bulk flow and parameters of the 
multipole components. The results of processing of the same 
subsamples in the non-relativ istic model are given in the paper 
jParnovskv & Parnowskill2010l) . 



2.3 Semirelativistic model of collective motion 

The obtained results appear to have problems - the values of S are 
200 times larger then their estimation from equation J 1 5 b - The same 
situation appeared when processing the previous sample. It is easy 
to see that the right-hand part of equation d 1 5b is dominated by the 
term containing tensor S due to the speed of light. 

For the homogeneous isotropic Universe the Weyl tensor and 
tensors S and Q vanish. They are connected with spatial inho- 
mogeneities of density distribution, e.g. attractors and voids. The 



left part of this relation is a sum of inputs of individual inhomo- 
geneities. For a spherically symmetric attractor with an excessive 
mass M at a distance u the spatial part of the tensor Cap^sV^V 5 
after reduction to eigenaxes receives the form: 



CafijsV V 



GM 














-1 



(17) 



This value falls cubically with distance, so the main 
input is provided by nearb y attractors. In the paper 
dParnovskv~ & Gavdamaka 2004) the input of the Great At- 
tractor, Perseum-Pisces superclaster, Shepley concentration and 
Virgo cluster were analyzed. The excess ive masses and distan ces 
to attractors were taken from the paper JMarinoni et al Jl 19980 . It 
was shown that the greatest input is provided by the Virgo cluster. 
Adding the inputs of all attractors we obtain an estimation of 
tensor S, which appears to be 200 times smaller than the calculated 
values. This is caused by the same distribution of the tensors Q 
and S over the celestial sphere. Due to measurement errors and 
deviations from the Tully-Fisher relation the tensor S "borrows" 
some of the value of tensor Q. Unfortunately, with the quality and 
quantity of available observational data we are unable to correctly 
separate the inputs of these tensors. Thus, taking into account the 
small values of tensor S , a so-called semirelativistic mo del was 
introduced in the paper jParnovskv & Gavdamakal |2004|) . which 
differs from the full relativistic model only by the dropped term 
with tensor S. So, it will be possible to use the relativistic model 
only when we get samples with significantly better quality and 
larger depth. For this reason, we switch to the semirelativistic 
model d2ir3ll4l [TTlll6t . The results of processing in semirelativistic 
model are also given in Table [T] The semirelativistic DQ- and 
D-models are exactly the same as their relativistic counterparts. 

However, in both the relativistic and semirelativistic models 
there is a serious problem. The calculated value of 7 appears to be 
way off the expected value, namely (—14.1 ± 2.9) ■ 10 -6 skm -1 
in the relativistic model and ( — 16.6 ± 2.6) • 10 -6 skm -1 in the 
semirelativistic model for Rmax = 10000 km s _1 . This value es- 
sentially differs from 70 = 3.98- 10~ 6 s km -1 dlOt calculated from 
cosmological parameters. As one can see from Table Q] the second 
term in equation J 14b is negligible and cannot be responsible for 
the discussed effect. Naturally, we do not question the values of 
the cosmological parameters and the reason should be sought else- 
where. 

In the next two section we consider the reasons, which could 
lead to the deviation of the calculated value of 7 from its value 
J 10b obtained from the cosmological parameters. These reasons can 
be caused either by the dependence of linear galaxy diameters on 
distance or by the influence of measurement errors. 



3 SHIFT OF 7 DUE TO DEPENDENCE OF LINEAR 
GALAXY DIAMETERS ON DISTANCE 

When determining the distances to the galaxies using the gener- 
alised Tully-Fisher relation in the "linear diameter - H I line width" 
version, we assume that galaxies with the same H I line width, mor- 
phological type, axial ratio and surface brightness index have the 
same linear diameter L. Let us consider the opposite case when the 
linear diameter weakly depends on the distance r according to the 
law 



L{r) = Lo + r _ 



(18) 
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Table 1. Parameters of the relativistic (R), semirelativistic (SR) and semirelativistic with fixed 7 (SR7) models of collective motion for the 
subsamples with Rmax = 10000 kms" 1 and Rmax = 8000 kms" 1 



R max = 10000 km s- 1 Rmax = 8000 km s" 1 







R 


SR 


SR7 


R 


SR 


SR7 


a, kms 1 




1130 


1134 


1157 


1008 


1018 


1033 


Ci 




-0.35 ±0.11 


-0.38 ± 0.11 


-0.39 ±0.11 


-0.62 ± 0.12 


-0.60 ±0.12 


-0.54 ±0.13 


C-2 




15.56 ± 1.26 


15.95 ± 1.25 


12.40 ± 1.18 


16.89 ± 1.43 


16.35 ± 1.42 


12.73 ± 1.31 


c 3 




1.96 ±0.19 


2.06 ±0.19 


l.bl±0.19 


2.55 ± 0.22 


2.59 ±0.22 


2.04 ±0.21 






9.23 ± 1.16 


9.45 ± 1.17 


8.06 ± 1.17 


8.84 ± 1.30 


9.73 ± 1.31 


8.45 ± 1.32 


c 5 




-1007 ± 103 


-1062 ± 103 


-566 ± 80 


-1095 ± 106 


-1120 ± 107 


-683 ± 81 


7, 10" 6 skm _ 


1 


-14.09 ±2.90 


-16.62 ± 2.55 


3.98 


-18.31 ± 3.85 


-19.59 ±3.59 


3.98 


D x , kms -1 




164.8 ± 93.8 


111.2 ± 91.7 


2.0 ±94.3 


154.8 ±99.3 


125.0 ±97.5 


42.7 ± 100.3 


D y , kms 




-117.4 ±95.2 


-191.8 ±90.1 


-108.6 ±93.2 


-67.8 ±95.0 


-100.2 ±93.4 


-113.0 ±96.6 


D z , kms" 1 




133.8 ± 76.5 


64.6 ± 72.0 


144.4 ± 73.9 


23.0 ± 76.6 


-44.4 ± 74.8 


32.4 ± 76.1 


D, kms" 1 




243 ± 79 


231 ± 80 


181 ± 74 


171 ± 90 


166 ± 89 


125 ± 87 


Z, b, deg 




325, 34 


300, 16 


271, 53 


336, 8 


321, -16 


291, 15 


gi, 10" 2 




4.90 ±4.28 


6.84 ± 1.39 


6.68 ± 1.65 


-3.33 ±5.11 


6.90 ± 1.41 


7.04 ± 1.67 


92, 10" 2 




5.76 ±4.66 


-1.99 ± 1.43 


-3.98 ± 1.70 


9.09 ±5.76 


-1.65 ± 1.49 


-1.89 ± 1.79 


93, 10- 2 




-3.11 ± 5.85 


-1.08 ± 1.82 


-0.47 ± 2.17 


-1.49 ± 7.39 


-1.10 ± 2.02 


-0.31 ± 2.41 


94, 10" 2 




-17.36 ±6.47 


2.22 ± 2.20 


3.79 ± 2.62 


-16.10 ± 7.78 


0.99 ± 2.17 


2.21 ± 2.58 


95, 10- 2 




-9.54 ±7.53 


0.37 ±2.40 


-0.65 ± 2.86 


-35.03 ± 9.08 


0.36 ± 2.41 


-0.47 ±2.87 


01 , 10 — 6 s km" 


-1 


3.18 ± 1.21 


2.35 ± 1.14 


3.40 ± 1.62 


3.79 ± 1.62 


3.34 ± 1.58 


4.55 ± 2.23 


02 , 10" 6 s km" 


-1 


-0.55 ± 1.52 


1.35 ± 1.32 


2.01 ± 1.87 


-4.05 ± 1.69 


-3.39 ± 1.61 


-5.08 ±2.26 


03, 10 — 6 skm" 


-1 


5.05 ± 1.45 


5.13 ± 1.37 


7.77 ± 1.95 


5.69 ± 1.85 


5.78 ± 1.82 


7.32 ± 2.57 


04, 10~ 6 skm" 


-1 


-4.31 ± 1.84 


-3.64 ± 1.70 


-4.82 ±2.41 


-3.75 ± 2.12 


-3.40 ± 2.09 


-5.46 ±2.96 


05, 10" 6 skm 


-1 


3.58 ± 1.60 


3.71 ± 1.49 


4.39 ±2.12 


3.44 ± 1.91 


3.82 ± 1.88 


4.28 ±2.65 


06, 10 -6 skm" 


-1 


-0.72 ± 1.69 


-1.86 ± 1.56 


-0.32 ± 2.21 


-3.20 ±2.04 


-3.56 ± 1.98 


-3.60 ±2.80 


07, 10 s km 


-1 


in n I r 01 

19.13 ± 5.81 


1 O 1 C 1 rr AC 

18. lb ± 5.4b 


21.03 ± 7.73 


on on 1 t tzr\ 

20.30 ± 7.59 


O A 1 £1 1 T A A 

24.16 ± 7.44 


31.06 ± 10.47 


pi, 10 -6 skm" 


-1 


-2.28 ± 1.74 


0.29 ± 1.48 


-2.49 ±2.09 


2.29 ±2.13 


4.66 ± 2.01 


3.11 ± 2.82 


P2, 10 -6 skm" 


-1 


-0.74 ± 1.94 


0.75 ± 1.79 


3.99 ±2.57 


-0.82 ± 2.65 


0.73 ± 2.55 


4.12 ± 3.67 


P3, 10 -6 s km" 


-1 


-1.75 ± 2.70 


1.41 ± 2.22 


-5.15 ±3.06 


-4.28 ±2.80 


-2.69 ± 2.65 


-5.95 ± 3.80 


si, 10~ 6 skm 


-1 


1.98 ± 5.55 






14.97 ± 7.41 






s 2 , 10" 6 skm" 


-1 


-10.88 ±5.82 






-16.09 ±8.28 






s 3 , 10" 6 skm" 


-1 


2.84 ±7.05 






1.29 ± 1.08 






s 4 , 10 -6 skm" 


-1 


28.01 ±8.53 






26.61 ± 1.15 






sb, 10 -6 skm" 


-1 


13.94 ±9.80 






54.83 ± 1.35 







Here Lo is the linear diameter of nearby galaxies according to 
Tully-Fisher relation. On the other hand, if we express the angu- 
lar diameters a in radians, the linear diameters will be given by 



L(r) = ar. 



(19) 



This yields a problem: we determine the distances to the galaxies 
using the Tully-Fisher relation in an assumption that the linear di- 
ameters of the galaxies are equal to Lo rather than L(r). Thus, the 
apparent distance p will be equal to 



L 

P = — • 

a 



(20) 



Combining the latter two formulae we can express the true distance 
r through the apparent distance p: 

Let us calculate the shift of 7 due to this effect. Radial veloci- 
ties of galaxies according to equation are given by 



V = Hr + 70 (Hrf, 



(22) 



where 70 is the true value of 7. However, instead of this formula 
we use the following expression: 

\2 



V = Hp + y(HpY 
It is trivial to find that 



(23) 



A7 = 7 — 70 = 



1 dL 
HLn ~~dr 



(24) 



There are several effects leading to the dependence l !24t . Let 
us consider them one by one. 



3.1 Effect of galaxy evolution 

Let us assume that galaxies evolve with time changing their linear 
diameters with a characteristic rate L = dL/dt. Since we observe 
more distant galaxies at earlier stages of development, we can write 



(25) 



dL _ _L _ J_L 

dr ~ c ' 1 ~~ HcV 

This effect is described, in particular, bv lWeinberjd ( ll972h . 

The question is, whether it alone can explain the observed 
value of A7 ~ 2.2 • 10 -5 s km -1 ? A simple estimation shows that 
would this be the case, the galaxies should shrink with a typical 
rate of 4.7 ■ 10~ 10 yr" 1 . This value contradicts to our knowledge 
of galaxy evolution. In particular, such large shrinking rate would 
lead to decrease of galaxies by about 12 per cent in one orbital pe- 
riod of the Sun in the Milky Way. Thus, this effect cannot make an 
essential contribution to the observed shift of 7. 
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3.2 Effect of selection by angular diameters 

The Tully-Fisher relation is statistical. In fact, galaxies can be 
smaller or larger than given by it. For galaxies with angular di- 
ameters close to the threshold at = 0'6 there is a selection. Large 
galaxies will enter RFGC and the small ones will have too small an- 
gular diameters and will be rejected. Thus, the average linear diam- 
eters of RFGC galaxies will increase with distance. This is usually 
referred to as Malmquist bias. According to equation d24t , these ef- 
fects will lead to a positive shift of A7. This shift has the opposite 
sign to the observed one and, therefore, cannot be its cause. Nev- 
ertheless, this effect cannot be totally neglected and further we will 
estimate its value using Monte Carlo simulations. 



3.3 Effect of cosmological decrease of observed surface 
brightness of galaxies 

Consideration of space-time curvature even in the simplest homo- 
geneous isotropic cosmological models leads to a number of ef- 
fects. One of these effects yields the decrease of observed surface 
brightness of galaxies. It is quite evident that the surface brightness 
is proportional to (r/D) 2 , where D is the photometric distance. In 
flat space-time there is no difference between r and D and the sur- 
face brightness does not depend on the distance. In curved space- 
time at small z we can use the ex pressions for r and D from the 
book dZeldovich & Novikovll 19831) , which yield I = I N r(1 - 4z). 
Here Inr is the surface brightness in the flat space-time and z 
is th e redshift. This expr ession is a low-z limit of the Tolman ef- 
fect jTolmanlll93dll934l) . which describes the decrease of surface 
brightness by a factor of (1 + z) 4 . Switching to the brightness 
/i = — 2.51g7 measured in mag/D" we obtain the shift 



A/X = p - UNR = J§Q R / C 

= 1.45 • 10 _d i?mag/D"skm- 1 . 



(26) 



For a galaxy at the distance lOO/i 1 Mpc this gives Afi = 
0.145 mag/D". Earlier such effects were considered, e.g. by 
ISandage&Lubrnl ( l200ll) . 

If galaxies had sharp edges where the surface brightness in- 
stantly vanishes, this effect would be of no interest for us. How- 
ever, for real galaxies the surface brightness gradually falls to the 
edges according to the de Vaucouleurs law 1(1) = Joe - '" , /x(Z) = 
Ho + 1.0857 Z/Zo ddeVaucouleurslI 19591) . Here lo and po are the 
values of I and /1 in the centre of the galaxy, I is the radial distance 
from the centre of galaxy and io is a characteristic radial scale. 
For spiral galaxies the scale Zo weakly depends on the morpho- 
logical type. For bright, comparable to the M ilky Way, galaxies it 
ranges from 1 to 10 kpc dReshetnikoy l l200lh . A more accurate es- 
timation was obtained by Fathi et alT i 2010h using the sample of 
more than 30000 galaxies. The average value of Zo appeared to be 
equal to 3.8 ±2.1 kpc. For smaller galaxies with the total stellar 
mass 10 9 - 10 10 Mq they estimated it as 1.5 ± 0.7 kpc and for 
larger galaxies with the total stellar mass 10 11 - 10 12 Mq - as 
5.7 ± 1.9 kpc. For the later types of galaxies, which are predom- 
inant in RF GC, the value o f Zo does not typically exceed 2 kpc, 
according to lFreemanl Jl970l) . In his sample of 36 galaxies the max- 
imum value of Zo was 6.1 kpc, and the second largest was 4.5 kpc. 

The linear diameters of RFGC galaxies are determined at 
the isophotal level /1 = 25 mag/D". Due to the abovementioned 
effect for distant galaxies this boundary is shifted with respect 
to its position in flat space-time. For instance, for a galaxy at 
r — lOO/i" 1 Mpc the region corresponding to the isophotal level 



25 mag/D " in the flat space-time would appear at the isophotal 
level 25.145 mag/D" in the real Universe. 

For this reason, the isophotal boundaries will be shifted to- 
wards the centre of the galaxy by SI = 1.33 • 10 _5 Zo^skm _1 . 
Naturally, the apparent decrease of isophotal diameter 5L will be 
twice this value. Using this distance-dependent isophotal diameter 
instead of constant linear diameter, we obtain the so-called isopho- 
tal distance instead of the angular diameter distance. 

The expression is derived for angular distance only. In re- 
ality, however, we deal with isophotal distances with a slightly dif- 
ferent value of 7. From the equations d26t and L(r) — Lo — SL 
we get ^7 = —2.66 • 10~ 5 Zo/I/o skm" . Let us estimate this 
value. As a typical galaxy diameter we take the value for the 
Milky Way: L = 30 kpc. The distribution of lg Lp for RFGC 
galaxi es, where Lo is expressed in kpc, is given by iKudrva et al.l 
(1997). The maximum of this distribution corresponds to the inter- 
val from 1.3 to 1.4, ehich corresponds to Lo ~ 22 kpc. The mean 
value of Lq for RFGC galaxies should be taken slightly larger be- 
cause (L > > 10 <lgLo) . With Lo = 30 kpc and l = 3 + 5 kpc 
we get A7 = (—2.7 -. — 4.4) • 10~ 6 skm _1 , which constitutes 
from 12 to 20 per cent of observed shift. Even the extreme es- 
timation with Lo = 22 kpc, Zo = 10 kpc can explain only 50 
per cent of the observed shift. For edge-on spiral galaxies the 
surface photometry was performed in a series of articles by van 
der Kr uit and Searle. The data are assembled in Table 3 of the 
paper ( Ivan der Kruit & Searldfl982l) . For 7 galaxies, 4 of which 
enter the RFGC (NGC4244=RFGC2245; NGC5907=RFGC2946; 
NGC4565=RFGC2335; NGC5023=RFGC2495), the ratio l /L 
ranges from 0. 10 to 0. 15, which corresponds to A7 = (—4 -. — 6) ■ 
lO^skm -1 . 

This effect can be responsible for the observed shift, but only 
partially. Thus, we still need to find the reason behind the main 
share of the shift. 



4 SHIFT OF 7 DUE TO THE INFLUENCE OF 
MEASUREMENT ERRORS 

Let us show that the observed shift can be explained with purely 
statistical effects due to measurement errors of H I line widths and 
angular diameters. 



4.1 Estimating the impact of measurement errors: a simple 
case 

Before trying to address this problem at its full extent, let us con- 
sider a simple case when an analytical solution can be provided. Let 
us start off from introducing a toy model y — Ax + Bx 2 , where 
x = W I a is the main term of the Tully-Fisher relation l |16l >, y = V, 
A = Ci, and B = 70C 2 . For generality we consider not only the 
value of 70 given by equation d 1 Ob but any fixed value. This model 
corresponds to an isotropic Hubble expansion with cosmological 
acceleration. 

Observational data provide us with a set of N points charac- 
terized by values Xi and yi . It is important to realize how measure- 
ment errors and deviations from the Tully-Fisher relation distort the 
dataset. The errors in velocity measurements and deviations from 
Tully-Fisher relation yield errors in y. The values of A and B deter- 
mined by the least square method from the dataset with such errors 
will have normal distribution of errors without shift. At the same 
time, errors in measurements of W or a yield errors in x. This case 
is similar to Malmquist bias. Due to this effect, the perceived values 
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of A and B given by the least square method will have a systematic 
error. 

Let us consider the following case: we have N values of n 
distributed uniformly over the interval [0, 1] withastep (jV — 1) . 
The values of y are calculated with A = Ao = 1 and B = Bo = 
7o. In each of N points the Xi is shifted by <r£j, where & is a nor- 
mally distributed quantity with zero mean and unit variance. The 
different values of f , are not correlated with each other. The values 
Hi are calculated from the original nondisplaced values of x%. The 
values of A and B given by the least square method have the form 



A = 



B 



JV JV 

E ViXi E xf 

8=1 i=l 



JV 2^3 

E E ^ 

8=1 8=1 



JV JV 

E E x? 

8=1 8=1 

JV 

„2„. „2 



E 



»=i 



JV JV JV JV 

E ^iVi E - E a=i E 

i— 1 i — 1 i — 1 i—1 



N 

E* 



JV 

2 E*f 

8=1 



JV 

E A 

8=1 



jV - 1 



' , > , i On 2 

Sj = + cr£i, yi = Xi + 70 (Xj ) . 



(27) 



(28) 



(29) 



It is not very difficult to calculate the mean values of A and B over 
£ using the following expressions: 



A = A x +70^2, 
B = Bi +70B2. 
Here we designated 



Ai = 



1 + 40a 2 + 60cr 4 



1 + 28a 2 + 180a 4 + 720<r 6 ' 



A 2 



28cr 2 + 60a 4 



1 + 28a 2 + 180a 4 + 720a 6 ' 



Si = 



S2 



-20a 2 + 120a 4 



1 + 28a 2 + 180a 4 + 720a B 
1 - 43/3a 2 + 80a 4 



(30) 
(3D 

(32) 
(33) 
(34) 
(35) 



1 + 28a 2 + 180a 4 + 720a 6 ' 

Note that the formulae l|32][33][34][35) are precise up to 0{1/N). 
To verify these formulae we calculated A and B for this toy model 
using 10000 Monte Carlo simulations. The mean values perfectly 
fitted the given formulae (see Figure QJ. 

If we calculate 7 using the formula 7 = B /A 2 , we will obtain 
instead of true value 70 a value 7, plotted on Figure[2]against a for 
different 70. At a = 0, i.e. when there are no errors, we obtain 
7 = 70, but at small a we obtain 7 < 70. This is the impact of the 
measurement errors we demonstrate. 




Figure 1. Dependence of the coefficients of equations i3Q) and 13 U on the 
noise level 




Figure 2. Dependence of least squares estimation of 7 on the noise level 
for different values of 70 . The values of 70 equal to values of 7 at <r = 



4.2 Estimating the impact of measurement errors: a real case 

The actual case is much more complicated. The problem can not 
be reduced to the one-dimensional case, since the galaxy's veloc- 
ity depends on its position on the celestial sphere. To estimate the 
distance we use all the terms in the relation d!6t . The errors in de- 
termination of distances are non-Gaussian. They are due to the er- 
rors in angular diameters and H I line widths and deviations from 
the Tully-Fisher relation. Thes e errors were analyzed in the paper 
dParnovskv & Parnowski 200§|). Here we will briefly mention the 
main points of the routine used. These errors can be described by 
four parameters: 



V i = VS°>(l + s v ^), (36) 
W> = W t (0) {l + s w &), (37) 
a l = af ) (l + s^ 3 ) + A a C4. (38) 

Here fi, £2, £3, and £4 are the four independent non-correlating 
random values. They are distributed according to Gauss law with 
zero mean and unit variance. 

Let us describe what errors correspond to each type of noise. 
The noise l |36t describes the deviations from the Tully-Fisher rela- 
tionship. It does not include velocity measurement error, because 
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radial velocities are well determined, as well as the directions to- 
wards galaxies. Such a type of noise provides a conventional log- 
normal distribution of the velocity deviation. If this noise is used 
alone, the maximal value of sv can be assumed about 0.2, which 
corresponds to 20 per cent uncertainty of distance estimated by the 
Tully-Fisher relation. 

The noise ( 1371 ) corresponds to H I line width measurement er- 
rors. Different methods of calculating Hi half-width of the same 
profile can differ by up to 10 per cent. Thus, if this noise is used 
alone, the maximal value of sw can be set to 0.15 with some toler- 
ance. 

The noise {38} describes angular diameter measurement er- 
rors. The value s a describes a relative error, caused by variations 
of exposition, curvature of galaxies etc. The diameter measure- 
ments have an error of about 5 ~ 10 per cent (Karachentseva, pri- 
vate communication). The value A a describes an absolute error 
of measurement. This value is important for the smallest galax- 
ies. Since these diameters were measured in tenths of millimetre, 
which corresponded to 0.11 arcmin, the value of A a can be esti- 
mated between 0.05 and 0.1. Note that if we use non-zero values of 
s a and A a simultaneously, the maximum estimation of A a should 
be somewhat reduced to avoid overestimating the noise added to 
angular diameters for the smallest galaxies. 

Naturally, such a difficult problem of determining the shift of 
7 can not be solved analytically. We use Monte Carlo simulations to 
resolve it numerically. Let us describe the details of this procedure. 

At first, we use a subsample with Rmax = 10000 km s _1 , 
which contains iV = 1459 galaxies. For each of the galaxies we 
substitute the measured radial velocity with the radial velocity cal- 
culated using the formulae i2l[3ll4l [TTl|16t . We use the values given 
by equations {8} and {9}. 

After that we add noise to our model and for each of 10000 
simulation we calculate the coefficients in the same way as we treat 
actual data, i.e. using the semirelativistic model. Thus, for each re- 
alisation we obtain a complete set of the coefficients including 7. 
For the obtained values of 7 we calculate the mean and the stan- 
dard deviation. Taking into account that the distribution of 7 is non- 
Gaussian, its quantiles differ from the ones calculated from the nor- 
mal distribution. In this article we give errors corresponding to the 
99 per cent confidence level. Application of Monte Carlo method 
allows to do this in a straightforward way. From 10000 values of 
7 for different realizations we find the 50th largest and smallest 
values. They give us the boundaries of the 99 per cent confidence 
interval. 

In principle, we can apply this procedure to any mock cata- 
logue. However, we should use a catalogue that has the same spatial 
distribution as well as distribution of morphological types, surface 
brightness index etc. as the main sample. The best mock catalogue 
is thus the sample itself. We use as a result the sample of actually 
measured parameters, namely angular diameter, surface brightness 
index, Hubble type, Hi line width, ratio of angular diameters in 
red and blue imprints, and celestial coordinates. The radial velocity 
is calculated from the formulae {7] [T6t with coefficients obtained 
from the real sample. Hwever, we use a more refined procedure to 
improve reliability. 

At first we calculated the coefficients for the semirelativistic 
model using real data while fixing the value of 7 at 70 ( |10t . The 
corresponding coefficients are also given in TableQ] Then we used 
this set of coefficients to calculate mock radial velocities for the 
galaxies. Such mock radial velocities are closer to the actual radial 
velocities than for any other models with fixed 7. 



The values obtained using the Monte Carlo method for dif- 
ferent parameters of the noise are given in Table [2] The top part 
of Table [2] illustrates the impact of each individual type of noise. 
One can see that the 7 value acts in the same way as in the simple 
case considered in the previous section. The noise sv, correspond- 
ing to the deviations from the statistical Tully-Fisher relation has 
little or no effect on the value of 7. All other types of noise lead to 
a drastic reduction of 7, especially the noise A a . The bottom part 
of Table[2]contains the results obtained with realistic noise param- 
eters. When choosing parameters we used as a control parameter 
the standard deviation a for the noised sample, comparing it to the 
corresponding value for real data, given in Table [TJ This param- 
eter is convenient because it vanishes in the absence of the noise 
and grows when the noise increases. Thus, it can help avoiding un- 
dernoising and overnoising. Other constr aints used and details of 
the p rocedure are described in the paper {Parnovskv & Parnowskil 
2008). Naturally, we do not try to find a unique set of noise pa- 
rameters. The four noise values form a four-dimensional parameter 
space, the a constraint yields a three-dimensional hypersurface in 
it. Other constraints give us rough estimates of the boundaries of 
the volume of suitable values. 

One can see from Table [2] that for a set of parameters lying 
inside or close to those boundaries, the value of 7 is much less 
than the initial value 70. Comparing these values with the value 
7 = (—16.6 ± 2.6) ■ 10 -6 s km -1 obtained from the real data, we 
can select a range of suitable noise parameters. For many realistic 
noise parameters the value of 7 falls into the la area, and for most 
of them it falls into 99 per cent confidence area. In those cases 
when 7 misses the 99 per cent confidence area, the value of A a is 
unrealistically large. 

We applied the same routine to the subsample with R max = 
8000 km s^ 1 . The results for this subsample are also presented in 
Tabled One can see that there are sets of noise parameters, which 
provide suitable shift of 7 for both subsamples. Nevertheless, we 
should mention that the same parameters provide a much less suit- 
able shift of 7 for the subsample with Rmax = 6000kms _1 , 
where the observed value is 7 = ( — 11.0 ± 7.6) ■ 10 -6 skm -1 . 
This is due to large errors in determination of 7 for subsamples 
with low depth. 

An additional advantage of this method is that by slight mod- 
ification of the algorithm we can also estimate the influence of se- 
lection described in section l3~2l For this purpose one should add 
an additional condition when adding noise to angular diameters: if 
the noised angular diameter becomes less than the minimal angular 
diameter in the actual sample, the noise should be reapplied. This 
procedure reduces the shift of 7, but this effect is not very strong, 
e.g. for Rmax = 10000 kms" 1 it gives 7 = (-15.6 ± 2.4) • 
10 -6 s km~ 1 against 7 = (-16.8 ± 2.3) ■ 10~ 6 s km^ 1 . For other 
subsamples this effect has the same order. The best noise parame- 
ters for this case are given in the bottom part of Table|2] 



5 ERRORS IN DETERMINATION OF DISTANCE 
CAUSED BY SHIFT OF 7 

An important result of shift of 7 is that the distances calculated 
by the relativistic Tully-Fisher relation l l 1 6t become more than the 
correct ones. Really, since the combination R + 7J? 2 is fixed by 
redshift data, the decrease of 7 leads to increase of R. We made 
some Monte Carlo simulations which showed that this increase is 
about 18 per cent. Such large errors will yield too large errors in 
determination of peculiar velocities, defined in a standard man- 
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Table 2. Results of Monte Carlo simulations, 3<r CL stands for 99 per cent confidence interval 



sir 



10000 kms" 



Rn 



8000 kms -1 



cr, kms 1 



7, 10 -6 skm" 



<7> ± cr 7 



3cr CL 



cr, kms 



7, 10- 6 skm _1 



<7> ± cr 7 



3ct CL 



Without correction for selection 



0.05 


0.00 


0.00 


0.00 


282 ±7 


4.0 ± 1.2 


( 1-1 


7.2) 


245 ±6 


4.0 ± 1.5 


( -0.9 


10.1) 


0.10 


0.00 


0.00 


0.00 


563 ± 14 


4.1 ± 2.4 


( -1.8 


10.4) 


491 ± 13 


4.1 ± 3.0 


( -6.7 


16.6) 


0.15 


0.00 


0.00 


0.00 


845 ± 21 


4.2 ± 3.6 


( —4.1 


14.3) 


■7*)/' 1 -* r\ 

736 ± 19 


4.3 ± 4.4 


(—12.8 


22.1) 


0.20 


0.00 


0.00 


0.00 


1126 ± 28 


4.4 ± 4.8 


( —6.7 


18.5) 


981 ± 25 


4.6 ± 5.9 


(—15.8 


30.8) 


0.00 


0.05 


0.00 


0.00 


312 ± 7 


—0.6 ±1.1 


( —3.4 


2.5) 


283 ± 7 


—3.2 ± 1.3 


( —8.0 


2.4) 


0.00 


0.10 


0.00 


0.00 


592 ± 14 


—9.4 ± 1.5 


(—12.8 


—5.2) 


531 ± 13 


—15.7 ± 1.6 


(—20.5 


—8.6) 


0.00 


0.15 


0.00 


0.00 


831 ± 19 


— 16.5 ± 1.3 


(—19.4 


— 12.8) 


nop? 1 -1 

737 ± 17 


—24.5 ± 1.2 


(—28.2 


— 18.8) 


0.00 


0.20 


0.00 


0.00 


1031 ± 23 


on n i n 

— 20.9 ± 1.0 


(—23.3 


—17.9) 


907 ± 21 


on a inn 

— 29.4 ± 0.9 


(—32.6 


—25.2) 


0.00 


0.00 


0.05 


0.00 


281 ± 7 


-2.2 ± 1.0 


( -4.8 


0.5) 


247 ±6 


-4.4 ± 1.2 


( -9.9 


0.0) 


0.00 


0.00 


0.10 


0.00 


538 ± 13 


-14.4 ± 1.2 


(-17.2 


-11.3) 


468 ± 11 


-19.8 ± 1.3 


(-24.2 


-14.2) 


0.00 


0.00 


0.15 


0.00 


756 ± 17 


-23.5 ± 0.8 


(-25.4 


-21.3) 


648 ± 15 


-30.2 ± 0.8 


(-32.7 


-25.8) 


0.00 


0.00 


0.20 


0.00 


946 ± 23 


-27.3 ± 1.0 


(-29.1 


-23.3) 


802 ± 21 


-33.8 ± 1.0 


(-36.2 


-21.6) 


0.00 


0.00 


0.00 


0.05 


370 ± 11 


-8.2 ± 1.4 


(-12.0 


-4.6) 


307 ± 10 


-10.4 ± 1.6 


(-16.9 


-4.3) 


0.00 


0.00 


0.00 


0.10 


664 ± 18 


-21.9 ± 1.2 


(-24.3 


-16.7) 


548 ± 15 


-27.2 ± 1.2 


(-31.3 


-1.0) 


0.17 


0.08 


0.03 


0.06 


1155 ±28 


-17.0 ±2.2 


(-23.6 


-5.9) 


1004 ± 25 


-22.6 ±2.4 


(-30.5 


-11.4) 


0.17 


0.06 


0.04 


0.06 


1125 ±27 


-16.0 ±2.3 


(-23.3 


-3.9) 


976 ± 25 


-21.1 ± 2.7 


(-28.8 


-9.1) 


0.17 


0.07 


0.04 


0.06 


1142 ±28 


-16.8 ± 2.2 


(-23.7 


-6.0) 


992 ± 25 


-22.2 ± 2.6 


(-29.8 


-11.2) 


0.18 


0.04 


0.02 


0.08 


1182 ±29 


-18.8 ± 2.1 


(-24.9 


-8.8) 


1020 ± 26 


-23.6 ± 2.4 


(-31.5 


-13.0) 


0.18 


0.05 


0.02 


0.07 


1169 ±29 


-16.8 ± 2.3 


(-23.5 


-5.4) 


1011 ± 26 


-21.5 ± 2.7 


(-30.0 


-9.7) 


0.18 


0.05 


0.03 


0.07 


1174 ±29 


-17.2 ± 2.3 


(-24.8 


-6.3) 


1016 ±26 


-22.0 ±2.7 


(-30.0 


-10.3) 


0.18 


0.05 


0.04 


0.07 


1181 ± 29 


-17.7 ±2.2 


(-24.9 


-5.9) 


1022 ± 26 


-22.8 ± 2.6 


(-30.1 


-11.1) 


With correction for selection 


0.18 


0.05 


0.02 


0.07 


1168 ±29 


-15.6 ± 2.4 


(-23.9 


-4.9) 


1009 ± 26 


-20.5 ± 2.8 


(-29.0 


-12.1) 


0.18 


0.05 


0.03 


0.07 


1173 ±29 


-16.0 ±2.4 


(-24.1 


-5.2) 


1014 ± 26 


-21.0 ±2.7 


(-29.2 


-13.0) 


0.18 


0.05 


0.04 


0.07 


1180 ±29 


-16.6 ± 2.3 


(-24.0 


-5.8) 


1020 ± 26 


-21.7 ±2.6 


(-30.5 


-13.8) 



ner V pec = V — R. This is a great drawback of relativistic and 
semirelativistic models. Note that after switching to curved space- 
time, it is possible to introduce an alternative definition of a pe- 
culiar velocity V pec = V — V CO sm, where V CO sm is the velocity 
of cosmological expansion, which is defined in the low-z limit as 
Vcosm = R + jR 2 . Peculiar velocities defined in this way suffer 
much less from this effect. Indeed, the non-relativistic model im- 
plicitly uses this second definition of peculiar velocity. If we take 
a look at the generalized Tully-Fisher relation @ we will see that 
the term with C's is, in fact, quadratic in distance. Thus, the 18 per 
cent difference in distances leads to underestimation of the quadru- 
pole components by a factor of 1.18, and of the octopole one - by 
a factor of (1.18) 2 . The dipole component remains unaltered. 

To correctly determine distances in relativistic or semirela- 
tivistic models we have two options. The first one is to introduce 
a correction for this effect. The second, technically easier, is to 
eliminate the origin of this effect by fixing the value of 7 at 70. 
In this way we use the information about cosmological constants 
obtained by other more precise methods. Thus we naturally switch 
to the next stage - the semirelativistic model with fixed value of 
7. Its 23 free parameters are given in Table [T] We also performed 
Monte Carlo simulations with the same noise parameters in this 
model. The difference between actual and estimated distances ap- 
peared to be about 0.25 per cent. Naturally, the same procedure can 
be applied to relativistic DQ- and D-models. In the next section we 
consider the collective velocity field obtained in the framework of 
this model of galaxy motion. 



6 THE MULTIPOLE STRUCTURE OF THE VELOCITY 
FIELD 

In this section we analyze the multipole structure of the velocity 
field. Nevertheless, we start from comparing the coefficients d 
of the generalized Tully-Fisher r elation to that obtained earlier by 
IParnovskv & Gavdamakal (|2004|). They changed not very signifi- 
cantly; one should note the decrease of the coefficient correspond- 
ing to the morphological type of the galaxy and a slight increase of 
the coefficient corresponding to the blue diameter. However, these 
trends are also present in the non-relativistic model and are caused 
by slightly different statistical properties of the updated sample. It 
is also interesting that the difference of the main coefficients Ci for 
the semirelativistic models with fixed 7 and with free 7 appeared 
to be 22 per cent, which is consistent with the 18 per cent differ- 
ence in distances for these two models and is naturally caused by 
the same reasons. 

For each regressor we calculated not only the coefficient and 
its error but also its statistical significance according to Fisher test. 
For the semirelativistic model with fixed 7 for the subsample with 
R max = 10000 kms -1 the minimum value F = 11.5 corre- 
sponds to Ci, the maximum value F = 109.7 - to C2. These 
values should be compared to the values 3.8, 6.6, 7.9, 10.8 and 
12.1, which correspond to 95, 99, 99.5, 99.9 and 99.95 per cent 
confidence levels respectively. Thus, all the coefficients of the gen- 
eralised Tully-Fisher relation dl6t are statistically significant at the 
99.9 per cent confidence level. 
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Figure 3. Bulk motion apices in galactic coordinates (Mollweide pro- 
jection) for Umax = 10000 km s —1 . Crosses mark the apices of 
the bulk motion in the D-model surrounded by la, 2a and 3a con- 
fidence areas. Solid boundaries correspond to results in the relativistic 
D-model with fixed 7, and dashed ones — to th e results of the non- 
relativistic D-model (Parnovsky & Parnowski 2010). Numbers denote the 
results of ot her authors: 1 - jLynde n-Bell et a l.lll988l). 2 - Hudson et alj 
Il995l). 3 - <Lauer&Postmanlll994l), 4 - jParnovskv et all l200ll). 5 
jDekel et alj l999h . 6 - Ida Costa et alj2000l) . 7 - |Hui 
IDale et al.ll 19991). 9 - feudrva et alj|2003l) . 10 - fWa 
- iParnovskv & Tugavl2004l) 
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Now let us consider the dipole component of the velocity field. 
Its parameters including the galactic coordinates I, b of the apex 
for the DQO-model are given in Table [T] The norms of the dipolar 
component do not contradict the ACDM model. For the model with 
fixed 7 the module of the dipolar component drops to 180 kms -1 . 
However, the bulk motion is usually considered in the framework 
of the simplest dipole models when the only characteristics of the 
velocity field are the modulus and the apex of the dipole compo- 
nent. In our case of DQO-models the velocity field is more com- 
plex and we cannot attribute the bulk motion solely to the dipole 
component. For this reason, to compare our results to the results 
of other authors we also calculated the dipole component in the 
framework of the relativistic (the same as semirelativistic) D-model 
with fixed 7. It yields the bulk flow velocity of 314 kms -1 di- 
rected towards I = 322°, b — 27° (Centaurus). On Figure [3] we 
plotted the boundaries of la, 2a and 3a confidence areas of this 
apex for Rmax = 10000 kms -1 . For this purpose we projected 
the 8-dimensional ellipsoid of errors into the 3-dimensional space 
and then projected it on the celestial sphere. On the same figure 
we also plotted the boundar ies of the confidence areas of t he apex 
in non-relativistic D-model ( Parnovsk y & Parnowskill2010h as well 
as positions of apices obtained by different authors. The value of 
the bulk motion appears to be larger than for DQO-models. For the 
subsample with R ma x = 8000 kms" 1 it is equal to 285 kms -1 . 
We see t hat D-models provide a result, which is closer to that ob- 
tained bv lWatkins et alj fc009T) . but still consistent with the ACDM 
model. 

Let us now consider the quadrupole component of the velocity 
field. What is the physical s ense of the quadrupole component? As 
one can see from the paper jParnovskv et al.|[200lh . it can be natu- 
rally combined with the Hubble constant. As a result, we obtain the 
effective 'Hubble constant' depending on direction 



H(l,b) = H(l + Qikrunk) 



(39) 



Naturally, this effective 'Hubble constant' is caused by the large- 
scale collective motion on the sample scale. To estimate the value 
of its anisotropy we found the eigenvalues and eigenvectors of ten- 



sor Q. The three eigenvectors are orthogonal and the sum of three 
eigenvalues is equal to zero because Q is a traceless tensor. 

We found the eigenvalues and the eigenvectors of the ten- 
sor Q for two considered subsamples. For the subsample with 
Rmax = 10000 kms -1 the maximal eigenvalue 7.1 ± 1.7 per cent 
corresponds to an axis directed towards I — 98° , b = 79° (Canes 
Venatici) and the opposite direction (Phoenix). The minimal eigen- 
value —4.1 ± 1.7 per cent corresponds to an axis directed towards 
/ = 195° ,6 = 2° (Gemini) and the opposite direction (Sagittar- 
ius). The third eigenvalue —3 per cent corresponds to an axis di- 
rected towards I = 286°, 6=11° (Centaurus-Vela) and the oppo- 
site direction (Androm eda-Lacerta). Comparing these values to the 
non-relativistic model jParnovskv & P amowski 2010) one can see 
that both the eigenvalues and the directions of the axes changed in- 
significantly. Nevertheless, the two negative eigenvalues, which are 
close to each other, have the opposite order in these two models. In 
this sense, the positive axis notably stands out, for which the effec- 
tive 'Hubble constant' exceeds the mean value by 7 per cent. For 
the subsample with Rmax = 8000 kms -1 the ellipsoid is three- 
axial and essentially differs from the oblate spheroid. The maximal 
eigenvalue 7.1 ± 1.7 per cent corresponds to an axis directed to- 
wards I — 102° , b = 85° (Canes Venatici) and the opposite direc- 
tion (Sculptor). The minimal eigenvalue —5.3 ±1.6 per cent corre- 
sponds to an axis directed towards I — 266° , b — 5° (Vela) and the 
opposite direction (Cygnus). The third eigenvalue —2 per cent cor- 
responds to an axis directed towards I = 356°, 6=1° (Sagittarius- 
Scorpio) and the opposite direction (Auriga-Taurus). These values 
are very close to those given by the non-relativistic model. Note that 
the axes for both subsamples nearly coincide with the exception of 
reverse order of negative eigenvalues for Rmax = 10000 km s -1 . 

We also calculated the statistical significance of these eigen- 
values. For Rmax = 10000 kms -1 the maximal eigenvalue has 
F = 18.2, which means that it is non-zero at 99.95 per cent confi- 
dence level, and the minimal eigenvalue has F = 5.8, which means 
that it is non-zero at 97.5 per cent confidence level. The similar sit- 
uation holds for the subsample with Rmax = 8000 kms -1 with 
Fisher values being 18.3 and 10.7 respectively. Additionally, we 
calculated the total statistical significance of the quadrupole com- 
ponent. The value y~ qua with its 5 degrees of freedom appears to be 
non-zero at over 99.5 per cent confidence level according to F-test. 

In the same way we calculated the total statistical signifi- 
cance of the octopole component. The value y oct with its 10 de- 
grees of freedom appears to be non-zero at over 99.5 per cent 
confidence level according to F-test. The value P with its 3 de- 
grees of freedom appears to be non-zero at slightly less than 90 
per cent confidence level according to F-test. Unlike the quadru- 
pole component, the octopole one lacks easily interpretable char- 
acteristics like eigenvector apices. The radial velocity field for 
R = 8000 kms -1 and Rmax = 10000 kms -1 in the semirela- 
tivistic model, which includes the octopole component, appeared 
to be very si milar to that in non-relativistic c ase, depicted on Fig. 6 
in the article jParnovskv & P arnowski 201fj|) . The most prominent 
feature of both these velocity fields is a strong inbound flow coming 
from the direction opposite to the apex of the bulk flow. 

Thus we obtained that the velocity field in the semirelativistic 
model with fixed 7 is very similar to that in non-relativistic case. 
The difference between these velocity fields may become signif- 
icant when more precise and deep samples will be available. For 
the existing sample this similarity yields two conclusions. The first 
one is that this similarity justifies the form of the generalised Tully- 
Fisher relation for the non-relativistic model. In contrast to the rel- 
ativistic model, the non-relativistic model ® was introduced em- 
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pirically. It includes the term, quadratic in distance, which has no 
theoretical substantiations. The relativistic model considered here 
explains why such a term needs to be included and what its or- 
der is. The second one is that the external information about the 
cosmological deceleration parameter, which we introduced into the 
model by fixing 7, is consistent with the observed parameters of 
RFGC galaxies. 

The bulk flow velocity is the most sensitive characteristic to 
the selection of the model used. It can change as much as 1.5 times 
when the same data is processed with different models of motion 
(314kms _1 forD-model, 281 kms -1 forDQ-model, lSlkms^ 1 
for DQO-model, and 249 km s" 1 for DQO-model without vector 
P - all with fixed 7 for IOO/1 -1 Mpc). This yields two consider- 
ations. First, the bulk flow velocity is a vulnerable characteristic 
of collective m otion for deep samples. Second, some authors like 
IWvmanl bold) believe that the excessively large values of bulk flow 
velocities obtained by some authors are a sufficient reason for aban- 
doning the ACDM cosmology for more exotic theories like brane 
cosmologies. We, however, have a different opinion on this mat- 
ter, and consider that such results should be double-checked using 
different models of collective motion. 



7 CONCLUSION 

We applied the relativistic model of motion supplied with the gen- 
eralised Tully-Fisher relation J 1 6b to the sample of 1623 flat edge- 
on spiral galaxies from the RFGC catalogue. The analysis of re- 
sults prompted us to switch first to the semirelativistic model, and 
then to the semirelativistic model with fixed 7. The parameters of 
the collective motion obtained in the framework of this model ap- 
peared to be close to that obtained in the non-relativistic case. We 
analysed certain reasons behind the decrease of 7 in the semirel- 
ativistic model. Evolution of galactic diameters, selection effects, 
and difference between isophotal and angular diameter distances 
appeared to be inadequate to explain this effect. At the same time, 
measurement error in Hi line widths and angular diameters can 
easily provide such a decrease. This was illustrated in a toy model, 
which allows analytical consideration, and then in the full model 
using Monte Carlo simulations. The obtained bulk flow velocity is 
consistent with ACDM cosmology. 



Kudrya Yu.N., Alexandrov A.N., 2002, J. Phys. Studies, 6, 472 
(in Ukrainian) 

Kudrya Yu.N., Alexandrov A.N., 2004, Herald Kyiv Univ. As- 

tron., 39-40, 130 (in Ukrainian) 
Kudrya Yu.N. et al., 1997, Astron. Lett., 23, 11 
Kudrya Yu.N. et al., 2003, A&A, 407, 889 
Lauer T.R., Postman M., 1994, ApJ, 425, 418 
Lynden-Bell D. et al., 1988, ApJ, 326, 19 
Marinoni C. et al., 1998, ApJ, 505, 484 

Parnovsky S.L., Gaydamaka O.Z., 2004, Kinematics and Physics 

of Celestial Bodies, 20, 477 
Parnovsky S.L., Tugay A.V., 2004, Astron. Lett., 30, 357 
Parnovsky S.L., Parnowski A.S., 2008, AN, 329, 864 
Parnovsky S.L., Parnowski A.S., 2009, preprint 

larXiv:0911.31021 
Parnovsky S.L., Parnowski A.S., 2010, Ap&SS, 325, 163 

larXiv:0910.46401 
Parnovsky SL. et al., 2001, Astron. Lett., 27, 765 
Reshetnikov V.P. "Surface photometry of galaxies", 
http://www.astro.spbu.ru/staff/resh/Book/index.html (in Rus- 
sian) 

Sandage A., Lubin L.M., 2001, ApJ, 121, 2271 
TolmanR.C, 1930, Proc. Nat. Acad. Sci., 16, 5111 
Tolman R.C., 1934, "Relativity, Thermodynamics, and Cosmol- 
ogy", Oxford Univ. Press: Oxford, 467 
van der Kruit PC, Seaiie L., 1982, A&A, 1 10, 61 
de Vaucouleurs G., 1959, Handbuch der Physik, 53, 311 
de Vaucouleurs G. et al., 1991, "Third Reference Catalogue of 

Bright Galaxies", Springer: Berlin, Heidelberg, New York 
Watkins R. et al., 2009, MNRAS, 392, 743 
Weinberg S., 1972, "Gravitation and Cosmology: Principles and 
Applications of the General Theory of Relativity", Wiley: New 
York, London, Sydney, Toronto 
Wyman M., Khoury J., 2010, preprint larXiv: 1004.20461 
Zeldovich Ia.B., Novikov I.D., 1983, "Relativistic astrophysics", 
University of Chicago Press: Chicago 

This paper has been typeset from a TgX/ LTj-X file prepared by the 
author. 



REFERENCES 

da Costa L.N. et al., 2000, ApJ, 537, L81 

Dale D.A. et al., 1999, ApJ, 510, Lll 

Dekel A. et al., 1999, ApJ, 522, 1 

Fathi K. et al., 2010, preprint larXiv: 1004. 15071 

Fisher R.A., 1954, "Statistical methods for research workers", 

Oliver and Boyd: London 
Freeman K.C., 1970, ApJ, 160, 811 

Hudson D.J., 1964, "Statistics Lectures on Elementary Statistics 

and Probability", CERN: Geneva 
Hudson M.J. et al., 1995, MNRAS, 274, 305 
Hudson M.J. et al., 2004, MNRAS, 352, 61 
Karachentsev I.D. et al., 1993, AN, 314, 97 
Karachentsev I.D. et al., 1999, Bull. SAO, 47, 5 

(astro-ph/0305566| > 
Karachentsev I.D. et al., 2000, Astron. Rep., 44, 150 
Komatsu E. et al., 2010, preprint larXiv: 1001.45381 



© 2010 RAS, MNRAS 000,[jJfTo] 



